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This document presents a measurement of W Z production in 1.02 fb 1 of pp collision data at ^/s = 7 TeV 
collected by the ATLAS experiment at LHC in 2011. A total of 71 candidates with a background expectation 
of 10.5±0.8(stat)l|^(sys) events were observed for purely leptonically decaying bosons with electrons, muons 
and missing transverse energy in the final state. The total cross section has been determined to be = 
21.1^2 g(stat) ^'^(syst) ^Q g(lumi) pb, in agreement with the Standard Model expectation of 17.2^J'g pb. 
Limits on anomalous triple gauge boson couplings have been derived. 



1. Introduction 

WZ Production in the standard model (SM) contains triple gauge coupling (TGC) vertex, which is essential to 
test the high energy behavior of the SM and is also a probe to new physics. Furthermore, new particles decaying 
into W ± Z pairs are predicted in supersymmetric models with an extended Higgs sector (charged Higgs) as well 
as models with extra vector bosons (e.g. W) pQ. 

At the LHC, the dominant W ± Z production mechanism is from quark- antiquark and quark-gluon initial 
states at leading-order (LO) and at next-to-leading order (NLO), respectively [2J. Only the s-channel diagram 
has a triple gauge boson interaction vertex and is hence the only channel to contribute to potential anomalous 
coupling behaviour of gauge bosons. 

This note presents a measurement of the W Z production cross section with the ATLAS detector in LHC 
proton-proton collisions at ^fs = 7 TeV using 1.02 fb -1 of data collected in 2011. The analysis uses four 
channels with leptonic decays (W Z — > ivli) involving electrons and muons: eee, ee\x, \x\ie or nnfj, plus missing 
transverse energy, _E™ 1SS . The main sources of background to the leptonic W Z signal are ZZ, Z7, Z+jets, 
and top-quark events. The signal and background contributions are modelled with Monte Carlo simulation and 
with data-driven measurements. 

Section [2] briefly describes the ATLAS detector and the data sample analysed in this paper. Section [3] 
discusses the signal and background simulation samples used in this analysis. The definition and reconstruction 
of physical observable objects such as particles and jets are detailed in Section |4j followed by event selection in 
Section [5j Section [6] presents the signal acceptance and background estimate, and the systematic uncertainties 
on these estimates. The calculations of the cross section and limits on the anomalous TGCs (aTGCs) are given 
in Section 



2. The ATLAS Detector and Data Sample 

The ATLAS detector is a multipurpose particle physics apparatus including a precision tracking system as 
the innermost part of the detector, highly segmented electromagnetic and hadronic calorimeters, and a large 
muon spectrometer [3]. 

This study uses data collected between March and June 2011. Data periods flagged with data quality 
problems that affect the reconstructed objects used in this analysis are removed. After data quality cuts, the 
total integrated luminosity used in this analysis is 1.02 fb -1 . The preliminary luminosity uncertainty for 2011 
data is 3.7%. 

W Z candidate events with multi-lepton final states are selected online with single- muon or single-electron 
triggers requiring px of at least 18 (20) GeV for muons (electrons). The trigger efficiency for W ± Z events, 
which pass all selection criteria, is in the range of 96-99% depending on the final state considered. 



3. Simulated Data Samples 

The W ± Z production processes and subsequent pure leptonic decays are modelled by the MC@NLO [1] gener- 
ator, which incorporates the NLO QCD matrix elements into the parton shower by interfacing to the Herwig [5] 
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program. The parton density function (PDF) set CTEQ6.6 [3] is used and the underlying event is modelled 
with Jimmy j7l[8j. The W^Z production cross section from qq' annihilation is calculated to be 17.2 jlg'| pb [4]. 
Electroweak corrections have not been considered as they are not relevant at the current luminosity [SI 110) . 

Major backgrounds for W ± Z signal detection are jets associated with W or Z gauge bosons, diboson and 
top-quark events. MC@NLO is used to model the tt and single top-quark events; to model the W/Z+]ets and 
Drell-Yan backgrounds we use Alpgen [TT] for e and /i decays and Pythia 12J for t decays. Events with heavy 
flavour dijets are modelled with PythiaB [T3]. The diboson processes WW and ZZ are modelled with Herwig 
and W/Z + 'y with MadGraph [14] , Whenever LO event generators are used, the cross sections are corrected by 
using k-factors to NLO or NNLO (if available) matrix element calculations [1]. The simulated background event 
samples generally correspond to an integrated luminosity of 1-10 fb _1 . Systematic uncertainties for simulated 
signal and backgrounds come from the uncertainty on their theoretical cross section, which are about 5-7% for 
diboson processes and 8-9% for tt. 

All event samples are simulated with multiple interactions in bunch crossings, including in-time and out-of- 
time pile-up, and with an average of 8 collision vertices per event. The number of interactions is re-weighted 
according to the luminosity distribution and the average number of interactions per bunch crossing of the data 
set used in this analysis, which varies according to the run period. 



4. Object Reconstruction 

The main physics objects necessary to select W^Z events are electrons, muons, and E^f iss . To ensure that 
the lepton candidates originate from a primary vertex, the longitudinal impact parameter with respect to the 
primary vertex must be less than 10 mm. Combined muons are identified by matching tracks reconstructed in 
the muon spectrometer (MS) to tracks reconstructed in the inner detector (ID). Their momentum is calculated 
by combining information from the two tracks and correcting for the energy loss in the calorimeter. ID tracks 
that were tagged as a muon on the basis of matching with a segment in the MS are also included. Only muons 
with px > 15 GeV and \r]\ < 2.5 are considered. The muon momentum in simulated events is smeared to 
account for a small difference between data and simulation. Non-prompt muons from hadronic jets are rejected 
by selecting only isolated muons, requiring the scalar sum of the track within AR < 0.2 of the muon 1 to be 
less than 10% of the muon p T [T5] , 

Electrons are formed by matching clusters found in the electromagnetic calorimeters to tracks in the inner 
detector. Electron candidates must have > 15 GeV, where Et is calculated from the cluster energy and 
track direction. To select central electrons and to avoid the transition regions between the calorimeters, the 
electron cluster must be within the regions I77I < 1.37 or 1.52 < |?7| < 2.47. Electrons are required to pass 
an electron identification requirement based on shower shape. To ensure isolation, the sum of the calorimeter 
energy in a cone of AR = 0.3 around the electron candidate, not including the energy of the candidate itself, 
must be less than 4 GeV. This quantity is corrected for the additional energy deposited in the presence of 
pile-up. The electron energy in simulated events is smeared to account for a small difference between data and 
simulation. Electrons overlapping with selected muons within AR < 0.1 are removed. 

The E™ ss is calculated with reconstructed electrons within \rj\ < 2.47, muons within \r/\ < 2.7, and jets and 
calorimeter energy clusters outside of other reconstructed objects within \rj\ < 4.5. The clusters are calibrated 
as electromagnetic or hadronic energy according to cluster topology. A small correction avoids double-counting 
the energy deposited by muons in the calorimeters [16] . 



5. Event Selection 

At least one single electron or muon trigger is required to fire in order to select the event, as described in 
Section [2j At least one primary vertex, with at least three good tracks associated, is required to remove non- 
collision backgrounds and ensure good object reconstruction. Events with two leptons of the same flavour and 
opposite charge with an invariant mass within 10 GeV of the Z boson mass are selected. This reduces much of 
the background from multi-jet and top-quark production, and some diboson backgrounds. 



x Ai? is defined as AR = V( Af ?) 2 + (A</> 2 ). 
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Events are then required to have at least three reconstructed leptons originating from the same primary 
vertex, two leptons from a Z — > 11 decay and one additional third lepton. The third lepton must pass more 
stringent identification criteria than the leptons attributed to the Z boson and have px > 20 GeV. The additional 
identification criteria required for electrons are on the matched track quality, the ratio of the energy measured 
in the calorimeter to the momentum of the matched track, and the detection of transition radiation. Muons are 
required to be reconstructed as combined muons. 

The transverse mass 2 of the W boson, M^F, is formed from the E™ ms and the third lepton. Furthermore, 
the Mp is required to be greater than 20 GeV. These cuts suppress the remaining backgrounds from Z and 
diboson production. For muon-triggered events, at least one of the muons is required to have px > 20 GeV and 
to have fired the muon trigger to ensure that the trigger is well onto the efficiency plateau above the threshold 
of the primary single muon trigger of 18 GeV. For electron-triggered events the px requirement is 25 GeV and 
the primary single electron trigger has a threshold of 20 GeV. 



6. Signal Acceptance and Background Estimate 

The acceptance of each cut is shown in Table []] after all corrections have been applied to the simulated W^Z 
events decaying leptonically to electrons and muons in the final state. 

Table I: MC W ± Z — > ivtt signal acceptance after each cut, and the relative acceptance is shown in parentheses. 



Cut Sequence 


Acceptance [%] 


Muon or electron trigger 


78.9 (78.9) 


Primary Vertex 


78.7 (99.8) 


mu - mz\ < 10 GeV 


28.2 (35.8) 


Three leptons 


12.3 (43.7) 


E^ iBB > 25 GeV 


10.0 (81.2) 


> 20 GeV 


8.5 (84.9) 


Trigger Match 


8.4 (99.5) 



The dominant backgrounds are events with jets associated with Z bosons, diboson production, and top-quark 
events. Backgrounds are estimated from a data-driven method where possible and from simulation otherwise. 

The data driven method uses a sample of Z events with an additional loose lepton which is dominated by Z+jet 
events from data. This sample contains two reconstructed leptons passing all particle identification requirements 
and one "lepton-like" jet that fails to satisfy lepton quality (medium e) or isolation (/i) requirements. The two 
leptons which pass all cuts must also pass the Z reconstruction requirements. To select a control sample as 
close to the signal region as possible, all other event selection criteria, including the E!f lss and Mp cuts are 
required. The event yield is obtained by scaling each event in the resulting sample by the "fake- factor" /(px), 
i.e. the probability that a "lepton-like" jet satisfies the quality or isolation requirements. The fake-factor is 
determined from a sample of events containing a Z boson plus an extra lepton-like jet, with a low missing energy 
requirement of E™ ss < 25 GeV, and extrapolated to high values of E™ lss using simulated events. The validity 
of the extrapolation has been verified with dijet events from simulation and data. 

The background from W/Z + 7 events where the photon converts into an electron-positron pair is not taken 
into account by the data-driven estimation methods, and is instead calculated with simulation. All other 
backgrounds are estimated using simulation. 

For the electron and muon objects, the uncertainties associated with the reconstruction and identification 
efficiency, energy scale, energy smearing, and isolation are taken into account. The uncertainties are determined 
from comparisons between simulated events and data in control samples and are around 2% to 6% depending 
on the decay channel. The uncertainties on the objects that are used to calculate _E™ 1SS are used to calculate the 
systematic uncertainties on E™ SB . Uncertainties in the description of the pile-up conditions by the simulation 
are also considered. The total systematic uncertainty on the acceptance of the E™ lss and transverse mass cuts 
due to the description by the simulation is 1-2%. 



2 The transverse mass Mx is defined as M^, = 2ExiEj' v — 2pxePT, 
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7. Results 

The numbers of expected and observed events after applying all selection cuts are shown in Tabic [TTJ Sta- 
tistical uncertainties are given for all four trilepton channels. 71 W ± Z candidates are observed in data, with 
10.5±0.8(stat)^2'i( s y s ) expected background events. The expected signal events include the contribution from 
r lepton decays into electrons or muons, which accounts for 1.7 events. The backgrounds from W + W~ and 
multi-jet production were found to be negligible. Various kinematic variables and the W ± charge of W ± Z 
candidates are shown in Figure [TJ 



Table II: Summary of observed events and expected signal and background contributions in the four trilepton and 
combined channels. Statistical uncertainties are shown for the individual channels, and both statistical and systematic 
uncertainties are shown for the data-driven background estimation methods are used for W/Z+]ets for all decay channels. 
The last row shows the ratio of expected signal events over the background expectation. 



Final State 

Observed 

ZZ 

W/Z+jets 



eee + h T ee/j, + E T 



11 

0.34±0.07 
2.03±0.38 



e/i/i + E^' BB fiLifi + ET Bt 
9 22 29 

1.03±0.13 0.82±0.12 1.40±0.15 
0.64±0.18 2.03±0.38 0.44±0.15 



Top 0.26T0.10 

W/Z + 7 0.49T0.28 

Total Background 3.08T0.49 

Expected Signal 7.55±0.17 
Expected S/B 2.5 



0.31±0.09 

I. 98T0.24 

II. 27±0.20 

5.7 



0.41±0.12 
0.56T0.39 
3.82T0.56 
12.12±0.22 
3.2 



0.60T0.15 

2.44T0.21 
18.16±0.27 
7.4 



combined 
71 

3.55±0.24±0.17 
5.14±0.59±^ 

1.58±0.23±0.10 
1.05±0.48±0.08 
10.5±0.8l2;? 
49.1±0.4±3.02 
4.3 
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Figure 1: Distributions of various kinematic variables, number of leptons, jets, and the W /± charge of W ± Z candidates. 
The kinematic variables shown are: (a) invariant mass of Z, (b) transverse mass of W ± , (c) transverse mass oiW ± Z, (d) 
transverse momentum of the Z (pf)> ( e ) £-™ lss , (f) lF ± charge, The points represent observed event counts with statistical 
errors, whereas the stacked histograms are the predictions from simulation including the statistical and systematic 
uncertainty. The last bin is an overflow bin. 
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7.1. Cross Section Extraction 

In order to combine the different channels, a common phase space region is defined in which a fiducial 
cross section is extracted. The common phase space is defined as p^ e > 15 GeV, \ri tl ' e \ < 2.5, > 25 GeV, 
\Mu — Mz\ < 10 GeV, and > 20 GeV, to approximate the event selection used in this analysis. 

For a given channel W ± Z — > lull where £ € {e, /i}, we define the fiducial cross section to be 

Ajobs _ ivrbkg AT MC 

o*« x B(WZ -> Ml) = a%U ivU = zTLzZl * (1 " Nff ) W 

Here, N^f e and N^j/ e denote the number of observed and background events respectively, C is the integrated 
luminosity and Cw z^Zvii is the efficiency corrected to account for measured differences in trigger and recon- 
struction efficiencies between simulated and data samples and for the extrapolation to the fiducial volume. 
B(W Z —> ivii) is the branching ratio for a W ± to decay to tv and a Z to decay to It, N^ IC is the number 
of W ± Z events where at least one of the bosons decays tauonically and N^f g c is the number of W ± Z events 
decaying into any lepton flavour. Since the fiducial volume is defined by the leptonic kinematics, the calculated 
cross section must include the branching ratio. The contribution from r lepton decays, about 4%, is not included 
and explicitly removed in the fiducial cross section definition with a scaling factor derived from the fraction of 
expected events from r decays. 

In addition to the fiducial cross section, we also calculate a cross section in the total phase space in each 
channel. The total cross section is calculated as: 

tot _ i^ivii ~ ^iuti) x (1 — ) (o) 

wz Cx B(WZ -> ivit) x Awz-ttM x C wz ^tuU N% g c ' { ' 

where Awz^iuU is the ratio of the number of events within the fiducial phase space region to the total number 
of generated events, needed to extrapolate to the total cross section for each channel, calculated using the 
event generator information. This ratio Awz^ivii equals 0.389, 0.392, 0.392, and 0.389 for the jU/i/i, e/j/j, ee/x, 
and eee channels, respectively. The uncertainties on Awz^iuU due to the MC signal sample statistics are 
estimated to be about 2% (uncorrelated among channels) and are included in the systematic uncertainties. The 
uncertainty due to the parton distribution functions is calculated by combining the uncertainty evaluated from 
44 CTEQ6.6 error eigenvectors (1.5%), and the central value deviation from MSTW2008 NLO (0.6%). The 
combined systematic uncertainty with quadratic sum is 1.6%. 

In the actual calculation, loglikelihood method is used to fit the cross-section and add nuisance parameters to 
account for the systematic uncertainty on the number of expected signal and background events in each channel, 
and to find the most probable value of a (fiducial or total) the negative log-likelihood function is simultaneously 
minimized over a and all the nuisance parameters Xk- 

The final results for the combined fiducial cross section for the W Z bosons decaying directly into electrons 
and muons (excluding contributions from r lepton decays) and the combined total inclusive cross section mea- 
surements are shown below. The fractional signal contribution from tau decays has been estimated with the 
simulated signal sample and is removed from the fiducial cross section 



a 



f ld -118+18/^^+6 



= 118±i°(stat) ±£(syst) l|(lumi) fb 



WZ^rtvii — - L - LU -16 



t$ z = 21.ll|;J(stat) iil(syst) ±^(lumi) pb. 



7.2. Limits on Anomalous Triple Gauge Couplings 

The model-independent method of effective action and effective Lagrangians allows for less restricted theories 
compared to the SM to be contructed. Expressions for the most general effective Lagrangian for a TGC vertex 
with two charged and one neutral vector boson can be found in Refs. |17j or [18] . If only terms that separately 
conserve charge conjugation C and parity P are considered then the Lagrangian reduces to 



i-wwz 
9wwz 



gf{W\ v W il Z v - W^W^Z") 



■ K Z WlW u Z" L 



X 



wt w^Z vp 

2 "* PM v " 



(3) 



where gwwz = — e cot Q^, 9w is the weak mixing angle, X^ u = d^X l/ — d^X fl and gf , k z and A are the coupling 
parameters that will be constrained in this analysis. It should be noted that gf , k z and A are dimcnsionless 
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coupling parameters that multiply dimension-four and dimension-six operators respectively; higher dimension 
operators are neglected. To recover the SM result <?f = 1, n z = 1, A = 0. In W Z production gf and A are 
proportional to s, whereas k z is proportional to \f§. The cross section has a quadratic dependance on a given 
coupling. 

To avoid tree-level unitarity violation at high centre-of-mass energies, which arises when radiative corrections 
from this new effective Lagrangian are bigger than the tree-level contributions, the anomalous couplings must 
vanish as Vs — > oo. To achieve this an arbitrary form factor is introduced according to 

a{s) = (4) 

where s is the invariant mass of the vector boson pair, a stands for Ag z , Ak z or A, the deviations of the 
anomalous couplings from the Standard Model values, and ao is the value of the anomalous coupling at low 
energy. A dipole form factor was used and A, the 'scale of new physics', was chosen to be 3 TeV. 



7.2.1. Reweighting and Limit Setting Procedure 

In MG@NLO version 4.0 [T5] it is possible to generate W Z events at any non-SM phase space point. Each 
event gets a vector of ten weights, wq — Wg, associated to it which allow an event sample to be reweighted to 
another phase space point. The weight at this new point is given by 



wtot (x wq + 2Ag z w\ + 2An z u)2 + 2\ws 

+2Ag z An z w i + 2Ag z \w 5 + 2An z Xw e 

+ (Ag z ) 2 w 7 + (A K z ) 2 w & + (X) 2 w 9 . (5) 

To set limits on aTGCs a profile likelihood test with a frequentist limit approach [5D] was adopted to determine 
all possible cross sections consistent with the observed data at the 95% confidence level (C.L.), giving upper and 
lower confidence limits for the cross section. The reweighting procedure allows the cross section to be determined 
as a function of aTGCs. Values of the aTGCs which predict cross sections inside of the 95% confidence interval 
(C.I.) of the cross section determine the 95% C.I. of the anomalous couplings. The C.I. for each anomalous 
coupling was determined separately with the other couplings set to their SM value. It should be noted that 
fiducial cross section values were used in this limit setting procedure, hence changes in Awz-^ivU are already 
taken into account. The effect of changes in Cwz-+tvtt within the range of couplings is ~1% and was not taken 
into account. 



7.2.2. Expected and Observed Limits 

Assuming SM values, the expected experimental limits on aTGCs can be found in Table \UJ\ where the median 
anomalous coupling C.I. is the best estimate, and the errors represent the 68% C.L. band of the anomalous 
coupling limits. The observed limits at the 68% and the 95% C.I on the anomalous couplings Agf , An z , and A 



can be found in Table IV For 68% C.I. two disjunct intervals are found. The observed 95% C.I. limits include 
the SM expectations and a comparison of ATLAS and Tevatron results are shown in Figure [2j LEP results from 
W + W~ production can be found in Ref. [3T]. Significant improvement in these limits is expected with more 
integrated luminosity and refined extraction methods taking advantage of the differential spectrum of kinematic 
quantities. 

The influence of the scale A on the observed limits was estimated by setting A above the centre-of-mass energy 
and calculating the expected cross section at the observed limits. The increase of the cross section is between 
2% and 7% depending on the energy dependence of the coupling parameters. This will result in tighter limits 
on the coupling values. The size of this effect is estimated to be about 10% on A, and less on An and gf . 



8. Conclusion 

A measurement of the W ± Z production cross section with the ATLAS detector in LHC proton- proton colli- 
sions at y/s — 7 TeV has been performed using final states with electrons and muons. With an integrated lumi- 
nosity of 1.02 fb _1 a total of 71 candidates is observed with a background expectation of 10.5±0.8(stat)i 2 ; 1 (sys). 
The Standard Model expectation for the number of signal events is 49.1±0.4(stat)±3.02(sys). The fiducial and 
total cross sections were determined to be 
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Table III: Expected experimental limits at the 95% C.I. on the anomalous couplings Apf , Ak z , and A, assuming SM 
values. The error bars indicate the 68% C.L. of the expected limits including systematic uncertainties for each coupling 
at the 95% C.L. 



Anomalous Coupling 


Limits of the 95% C.L 


Ag? 


[-0.16^,0.24^] 


Ak z 


h0.7l° ; 2 2 ,0.9; ° 2 2 ] 


A 


[-0.14±g;g*,0.14^g : g*] 



Table IV: Observed limits at the 68% and 95% C.L on the anomalous couplings Apf , A« z , and A. 



Anomalous Coupling 


Limits of the 68% C.L 


Limits of the 95% C.L 


Afff 


[-0.17,-0.05] , [0.13,0.26] 


[-0.21,0.30] 


Ak z 


[-0.8, -0.2], [0.5, 1.0] 


[-0.9,1.2] 


A 


[-0.15,-0.06] , [0.06,0.15] 


[-0.18,0.18] 



°wz^h,u = H8li?(stat) ^(syst) t|(himi) fb and d% z = 21.ll^(stat) j£l(syst) tgllflumi) pb, respectively. 
The total cross section is in good agreement with the SM expectation g] of 17.2l ; 8 pb. Limits on the anomalous 
triple gauge couplings gf , k z and A are derived. 



X 



1 1 1 1 1 1 1 1 1 1 1 1 

ATLAS Preliminary 


i i i i | i i i i | i i 
ATLAS. JYdt = 1 .0 fb 1 
Vs=7TeV, A=3TeV — 

__ CDF, ^Ldt = 1 .9 fb' 1 

Vs=1 .96TeV, A=2TeV 

DO. ^Ldt = 4.1 fb' 

Vs=1.96TeV,A=2TeV — 






W ± Z^> Nil 
95% C.L. 

i i 1 i i i i 1 i i i i 


i i i i I i i i i I i i 



-1 -0.5 0.5 1 



Figure 2: Shown are aTGC limits from ATLAS and Tevatron experiments. CDF [22 and D0 [23] limits are for W ± Z 
production with a Pt(Z) shape fit; ATLAS limits are for a cross section fit. Luminosities, centre-of-mass energy and 
cut-off A for each experiment are shown and the limits are for 95% C.L. 
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